Toward an Ising Model of Cancer and Beyond

Theoretical and computational tools that can be used in the clinic to predict neoplastic progression and propose individualized optimal treatment strategies to control cancer growth is desired. To develop such a predictive model, one must account for the complex mechanisms involved in tumor growth. Here we review resarch work that we have done toward the development of an "Ising model" of cancer. The review begins with a description of a minimalist four-dimensional (three in space and one in time) cellular automaton (CA) model of cancer in which healthy cells transition between states (proliferative, hypoxic, and necrotic) according to simple local rules and their present states, which can viewed as a stripped-down Ising model of cancer. This model is applied to model the growth of glioblastoma multiforme, the most malignant of brain cancers. This is followed by a discussion of the extension of the model to study the effect on the tumor dynamics and geometry of a mutated subpopulation. A discussion of how tumor growth is affected by chemotherapeutic treatment is then described. How angiogenesis as well as the heterogeneous and confined environment in which a tumor grows is incorporated in the CA model is discussed. The characterization of the level of organization of the invasive network around a solid tumor using spanning trees is subsequently described. Then, we describe open problems and future promising avenues for future research, including the need to develop better molecular-based models that incorporate the true heterogeneous environment over wide range of length and time scales (via imaging data), cell motility, oncogenes, tumor suppressor genes and cell-cell communication. The need to bring to bear the powerful machinery of the theory of heterogeneous media to better understand the behavior of cancer in its microenvironment is presented.Comment: 55 pages, 21 figures and 3 tables. To appear in Physical Biology. Added reference

Mathematical and Statistical Techniques for Systems Medicine: The Wnt Signaling Pathway as a Case Study

The last decade has seen an explosion in models that describe phenomena in systems medicine. Such models are especially useful for studying signaling pathways, such as the Wnt pathway. In this chapter we use the Wnt pathway to showcase current mathematical and statistical techniques that enable modelers to gain insight into (models of) gene regulation, and generate testable predictions. We introduce a range of modeling frameworks, but focus on ordinary differential equation (ODE) models since they remain the most widely used approach in systems biology and medicine and continue to offer great potential. We present methods for the analysis of a single model, comprising applications of standard dynamical systems approaches such as nondimensionalization, steady state, asymptotic and sensitivity analysis, and more recent statistical and algebraic approaches to compare models with data. We present parameter estimation and model comparison techniques, focusing on Bayesian analysis and coplanarity via algebraic geometry. Our intention is that this (non exhaustive) review may serve as a useful starting point for the analysis of models in systems medicine.Comment: Submitted to 'Systems Medicine' as a book chapte

Two-scale Moving Boundary Dynamics of Cancer Invasion:Heterotypic Cell Populations Evolution in Heterogeneous ECM

This book contains a collection of original research articles and review articles that describe novel mathematical modeling techniques and the application of those techniques to models of cell motility in a variety of contexts. The aim is to highlight some of the recent mathematical work geared at understanding the coordination of intracellular processes involved in the movement of cells. This collection will benefit researchers interested in cell motility as well graduate students taking a topics course in this area.

Expression analysis of E-cadherin, Slug and GSK3β in invasive ductal carcinoma of breast

<p>Abstract</p> <p>Background</p> <p>Cancer progression is linked to a partially dedifferentiated epithelial cell phenotype. The signaling pathways Wnt, Hedgehog, TGF-β and Notch have been implicated in experimental and developmental epithelial mesenchymal transition (EMT). Recent findings from our laboratory confirm that active Wnt/β-catenin signaling is critically involved in invasive ductal carcinomas (IDCs) of breast.</p> <p>Methods</p> <p>In the current study, we analyzed the expression patterns and relationships between the key Wnt/β-catenin signaling components- E-cadherin, Slug and GSK3β in IDCs of breast.</p> <p>Results</p> <p>Of the 98 IDCs analyzed, 53 (54%) showed loss/or reduced membranous staining of E-cadherin in tumor cells. Nuclear accumulation of Slug was observed in 33 (34%) IDCs examined. Loss or reduced level of cytoplasmic GSK3β expression was observed in 52/98 (53%) cases; while 34/98 (35%) tumors showed nuclear accumulation of GSK3β. Statistical analysis revealed associations of nuclear Slug expression with loss of membranous E-cadherin (p = 0.001); nuclear β-catenin (p = 0.001), and cytoplasmic β-catenin (p = 0.005), suggesting Slug mediated E-cadherin suppression via the activation of Wnt/β-catenin signaling pathway in IDCs. Our study also demonstrated significant correlation between GSK3β nuclear localization and tumor grade (p = 0.02), suggesting its association with tumor progression.</p> <p>Conclusion</p> <p>The present study for the first time provided the clinical evidence in support of Wnt/β-catenin signaling upregulation in IDCs and key components of this pathway - E-cadherin, Slug and GSK3β with β-catenin in implementing EMT in these cells.</p

Cell-scale degradation of peritumoural extracellular matrix fibre network and its role within tissue-scale cancer invasion

Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. Occurring over many different temporal and spatial scales, the first stage of invasion is the secretion of matrix degrading enzymes (MDEs) by the cancer cells that consequently degrade the surrounding extracellular matrix (ECM). This process is vital for creating space in which the cancer cells can progress and it is driven by the activities of specific matrix metalloproteinases (MMPs). In this paper, we consider the key role of two MMPs by developing further the novel two-part multiscale model introduced in [33] to better relate at micro-scale the two micro-scale activities that were considered there, namely, the micro-dynamics concerning the continuous rearrangement of the naturally oriented ECM fibres within the bulk of the tumour and MDEs proteolytic micro-dynamics that take place in an appropriate cell-scale neighbourhood of the tumour boundary. Focussing primarily on the activities of the membrane-tethered MT1-MMP and the soluble MMP-2 with the fibrous ECM phase, in this work we investigate the MT1-MMP/MMP-2 cascade and its overall effect on tumour progression. To that end, we will propose a new multiscale modelling framework by considering the degradation of the ECM fibres not only to take place at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour

Development of a Three Dimensional Multiscale Computational Model of the Human Epidermis

Transforming Growth Factor (TGF-β1) is a member of the TGF-beta superfamily ligand-receptor network. and plays a crucial role in tissue regeneration. The extensive in vitro and in vivo experimental literature describing its actions nevertheless describe an apparent paradox in that during re-epithelialisation it acts as proliferation inhibitor for keratinocytes. The majority of biological models focus on certain aspects of TGF-β1 behaviour and no one model provides a comprehensive story of this regulatory factor's action. Accordingly our aim was to develop a computational model to act as a complementary approach to improve our understanding of TGF-β1. In our previous study, an agent-based model of keratinocyte colony formation in 2D culture was developed. In this study this model was extensively developed into a three dimensional multiscale model of the human epidermis which is comprised of three interacting and integrated layers: (1) an agent-based model which captures the biological rules governing the cells in the human epidermis at the cellular level and includes the rules for injury induced emergent behaviours, (2) a COmplex PAthway SImulator (COPASI) model which simulates the expression and signalling of TGF-β1 at the sub-cellular level and (3) a mechanical layer embodied by a numerical physical solver responsible for resolving the forces exerted between cells at the multi-cellular level. The integrated model was initially validated by using it to grow a piece of virtual epidermis in 3D and comparing the in virtuo simulations of keratinocyte behaviour and of TGF-β1 signalling with the extensive research literature describing this key regulatory protein. This research reinforces the idea that computational modelling can be an effective additional tool to aid our understanding of complex systems. In the accompanying paper the model is used to explore hypotheses of the functions of TGF-β1 at the cellular and subcellular level on different keratinocyte populations during epidermal wound healing

3D Multi-Cell Simulation of Tumor Growth and Angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors

Structured models of cell migration incorporating molecular binding processes

The dynamic interplay between collective cell movement and the various molecules involved in the accompanying cell signalling mechanisms plays a crucial role in many biological processes including normal tissue development and pathological scenarios such as wound healing and cancer. Information about the various structures embedded within these processes allows a detailed exploration of the binding of molecular species to cell-surface receptors within the evolving cell population. In this paper we establish a general spatio-temporal-structural framework that enables the description of molecular binding to cell membranes coupled with the cell population dynamics. We first provide a general theoretical description for this approach and then illustrate it with two examples arising from cancer invasion